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We analyze in detail the peculiar velocity field traced by 56 clusters within 120 h 1 Mpc 
in the "Streaming Motions of Abell Clusters" (SMAC) sample. The bulk flow of the SMAC 
sample is 687 ± 203 km s _ \ toward I = 260° ± 13°, b = 0° ± 11°. We discuss possible 
[ systematic errors and show that no systematic effect is larger than half of the random error. 

■<^j- ■ The flow does not drop off significantly with depth, which suggests that it is generated by 

f^S structures on large scales. In particular, a Great Attractor as originally proposed by Lynden- 

Bell et al. cannot be responsible for the SMAC bulk flow. The SMAC data suggest infall into 
■ an attractor at the location of the Shapley Concentration, but the detection is marginal (at the 

90% confidence level). We find that distant attractors in addition to the Shapley Concentration 
are required to explain the SMAC bulk flow. A comparison with peculiar velocities predicted 
from the IRAS PSCz redshift survey shows good agreement with a best fit value of (3\ = 
^m°' 6 /&i = 0.39 ± 0.17. However, the PSCz density field is not sufficient to acount for 
all of the SMAC bulk motion. We also detect, at the 98% confidence level, a residual bulk 
flow of 372 ± 127 km s _1 toward I = 273°, b = 6° which must be generated by sources not 
included in the PSCz catalogue, that is, either beyond 200 /i _1 Mpc, in the Zone of Avoidance 
or in superclusters undersampled by IRAS. Finally, we compare the SMAC bulk flow with 
other recent measurements. We argue that, at depths ranging from 60 to 120 /i _1 Mpc, flows 
of order 600 km s _1 are excluded by multiple data sets. However, convergence to the CMB 
& ' frame by a depth of 60 h~ 1 Mpc is also excluded by multiple data sets. We suggest that a bulk 

flow of 225 km s _1 toward I = 300°, b = 10° at depths greater than 60 /i _1 Mpc is consistent 
with all peculiar velocity surveys, when allowance is made for sparse sampling effects. 
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ies: clusters: general - surveys - cosmological parameters - large-scale structure of Universe 



1 INTRODUCTION 

The source of the Local Group's (LG) peculiar velocity of 627 ± 
22 km s" 1 toward I = 276°, b = 30° <Kogut et alJll993l) with re- 
spect to the Cosmic Microwave Background (CMB) has been a 
puzzle since the detection of the CMB dipole ISmoot et alll977t) . 
In order to resolve this fundamental question, it is necessary to map 
the peculiar velocities of nearby galaxies. If a volume which is at 
rest with respect to the CMB can be identified, then the masses re- 
sponsible for the motion of the LG must be contained within that 
volume. 

An important milestone in p eculiar velocity surveys was the 
study of iLvnden-Bell et al] 1 1988), who used the D n — a distance 
indicator and claimed detection of a "Great Attractor" (GA) at a 
distance of approximately 45 /i _1 Mpc believed to be responsi- 



ble for most the LG's motion. Subsequent peculiar velocity sur- 
veys of the nearby Universe undertaken in the early 90s did not 
reveal the expected in fall signature on the far side of the GA 
I Mathewson et al. 1992) and hinted at a large coherence length 
for the flow <Willicldll99(ilCourteadll992h Analyses of redshif t 
surveys <DressleJll988t IStrauss et alJll992l: lHudsorj[l993l 1 1994 
did identify a significant excess of galaxies near the proposed GA, 
but with an abundance insufficient to generate the LG's motion. 
For example, Hudson <1994 concluded that 400 ± 45 km s -1 of 
the LG's motion arose from sources beyond 80 /i _1 Mpc. The lo- 
cation of t he peak of the GA in the Zone of Avoidance (ZOA), 
as given bv lKolatt et alj Il995t) . suggests the possibility that most 
of the GA's overdensity may be obscured. Searches in the ZOA 
have r evealed the presence of massive clusters such as A3627 
I Kraa n-Kortewee et alll99(j) but the integrated overdensity within 
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the ZOA appears still insufficient to generate substantial infall at 
the LG IStavelev-Smith et all2000h . 

lLauer & Postmarj ll994l hereafter LP) used the photomet- 
ric properties of brightest cluster galaxies as standard candles 
and claimed that a much larger region of space extending to 150 
/i _1 Mpc was moving at a velocity of 689 ± 178 km s _1 toward 
I = 343°, b = 52°. The LP result was puzzling, first because the 
amplitude of the flow on such a large s cale wa s higher than ex- 
pected in popular cosmological models IStrauss et alJll995h . and 
second because the direction of the LP flow was significantly dif- 
ferent from previous bulk flow measurements l Stra uss & Willicla 
Il995t) . 

Recent surveys of field galaxies to depths R < 60/i _1 Mpc 
generally suggest that the bulk flow within these nearby volumes 
is in the range from 100 - 300 km s -1 in the direction I ~ 300°, 
b fts 10° Jwillick et alll99ilGiovanelli et alll998tlCourteau et all 
l2000Hda Costa et all200CtlTonrv et all200d) . Clusters surveys on 
large (R > lOO/i -1 Mpc) scales have not supported the LP result, 
but, at face value, also app ear to suggest a wide range of values 
for the large- s cale motion (M ould et al|l993j:|H udson et alll999t 
IWillickll999tlDale et aJl999tlColless et all200ll) . This apparent 
disagreeme nt is due to neglect of spa rse sampling on the quoted 
flow errors JWafkins & Feldmarl l995). When the effects of sparse 
sampling are properly taken into account, the large-scale surveys 
are not in conflict (Hudson 2003). For re cent reviews, th e reader is 
directed to lCourteau & Dekel J200ll) and lZaroubil 120021) . 

This paper is the fifth in a series based on the "Streaming Mo- 
tions of Abell Clusters" (SMAC) project. The aim of this project 
is to obtain distance estimates, via the Fundamental Plane (FP) 
method, for elliptical galaxies in clusters and to map the peculiar 
velocity field within 120 fe -1 Mp c of the Local Group. First results 
from this work were presented in Huds on et alj |1999>. who quoted 
a bulk flow of 630 ± 200 km s" 1 toward I = 260°, b = -1° from 
a preliminary analysis of the same sample studied here. 

Pr evious papers in this serie s have reported new data for this 
project. ISmifh et al J <200(ll200lL hereafter SMAC-I and SMAC-II 
respectively) presented spectroscopic and photometric data, respec- 
tively; these data were compared and combined with other data 
availabl e in the literature to obtain the final data set used for this 
project iHudson et alll200lL hereafter SMAC-III). Smith et al. (in 
preparation, hereafter SMAC-IV) reports FP distances for 56 clus- 
ters. 

In this paper, we analyze the peculiar velocity field in the Uni- 
verse as traced by the SMAC clusters. In Section|2| we summarize 
the SMAC peculiar velocity data. In Section[3] we model the pecu- 
liar velocity field as a simple bulk flow (the dipole moment of the 
velocity field). An important aspect of this paper is to evaluate the 
robustness of the results to systematic effects (Section f3.3> . In Sec- 
tion |4| we consider more complicated flow models, including toy 
models based on simple attractors. Section|5|compares the peculiar 
velocities of SMAC clusters to the predictions from the IRAS PSCz 
density field. Section[6]compares the flow of the SMAC sample to 
other results in the literature. 



2 DATA 

In SMAC-IV, we tabulated peculiar velocity data for 56 clus- 
ters spanning the whole sky, and extending to a distance of ~ 
120ft -1 Mpc. The peculiar velocities quoted there were based on a 
inverse fit to the FP, in which log a is regressed on log R c and (/i) c . 



The inverse fit has the advantage that it is insensitive to selection 
on the photometric parameters. 

Median peculiar velocity errors per cluster are ~ 575 km s~ . 
These errors are primarily due to intrinsic scatter in the Fun- 
damental Plane which contributes a fractional distance error of 
0.21/\/iV, where TV is the number of galaxies in the cluster. How- 
ev er, the errors also inc lude the effect of the 16% uncertainties in 
the Schleael et al. 1 1998, hereafter SFD) extinctions and uncertain- 
ties in the mean redshifts of the clusters. 

In SMAC-IV, we also quoted for each cluster the error contri- 
bution arising from uncertainties in matching velocity dispersions 
from disparate observing runs. An important aspect of this work is 
the treatment of these errors, which are not independent from clus- 
ter to cluster, but are correlated on the sky. To calculate their effect 
on the bulk flow, we bootstrap re-sample the "overlap" sample used 
to calculate the velocity dispersion matching corrections. We then 
generate new FP datasets using those bootstrapped corrections and 
use these bootstrap samples to calculate the covariance matrix of 
the cluster peculiar velocities. We find that these systematic errors 
do not dominate the error budget in the bulk flow. On the other 
hand, due to their coherent nature, neither can they be neglected 
(see Section U.3.1l below). 



3 THE BULK FLOW 
3.1 Flow Model 

One statistic of particular interest is the bulk flow, or the dipole 
moment of the peculiar velocity field. For an idealized densely- 
sampled survey, the bulk flow would reflect the gravitational pull of 
mass near to and beyond the survey limits. It is therefore sensitive 
to the distribution of mass on the largest scales. 

We fit the radial components of the cluster peculiar velocities 
with the flow model 

A TJ 

U(r) = V-f + —r (1) 
ti 

where the free parameters are V, the bulk flow vector, and AH/H, 
a perturbation to the assumed Hubble constant. The distances 
quoted in SMAC-IV are in units of km s _1 and have already been 
adjusted so that the best fit gives AH = 0, but leaving this param- 
eter free allows for the correct propagation of errors into the bulk 
flow V. We use Galactic coordinates as the basis for our Cartesian 
coordinates: z is toward the North Galactic Pole (6 = 90°), x is 
in the direction of the Galactic centre (I — 0°), and y is in the 
direction of rotation (I — 90°). 

As noted above, the errors in the peculiar velocities are not 
independent but rather are coupled through the velocity-dispersion 
system-matching corrections described above and in Paper III. We 
construct a 56 x 56 covariance matrix which consists of the variance 
in distance due to scatter in the FP plus the square of a "thermal" 
scatter <j t h = 250 km s _1 , where the latter term reflects the small- 
scale noise of individual clusters around the mean bulk flow. These 
errors are independent from cluster to cluster so this part of the 
covariance matrix is diagonal. This is then added to the system- 
matching covariance matrix to obtain the full covariance matrix, C 
. We then minimize 

X 2 = - U(n)] [C-% [u 3 ~ U(r 3 )] (2) 

where the C" 1 is the inverse covariance matrix. 



Streaming motions of galaxy clusters — V. The peculiar velocity field 3 



To assess the depth of a sparse sample such as SMAC, we 
define an error-weighted effective depth, 

_ Ei r V^i 

R* — , , 5- (3) 

where of are the peculiar velocity errors, i.e., the diagonal elements 
of C. For an idealized, densely-sampled survey filling a top-hat 
sphere of radius R with uniform distance errors for all objects, one 
would obtain R t = R/2. With our choice of error weighting, the 
SMAC sample has R, ~ 6300 km s , as expected for a survey 
extending to ~ 12000 km s -1 . 



3.2 Results 

For our preferred choice <r t h = 250 km s , the SMAC bulk flow 
is 687 ± 203 km s"\ toward I = 260° ± 13°, b = 0° ± 11°. In 
Galactic Cartesian coordinates, this is V x = —124 ± 172 km s _ , 
V y = -676 ± 190 km s" 1 , V z = -5 ± 131 kms~\ and in 
Supergalactic Cartesian this becomes Vsgx = —367 kms~ , 
Vsgy = 59 kms" 1 , Vsgz = —578 kms^ 1 , in the direction 
SGL = 171°, SGB = -57°. This flow is significantly different 
from zero at the 99.91% confide nce level (CL). In some papers 
lLauer & Postmanl 19941 including Hu dson et al J 19991) . the authors 
allow for "error-bias". This is a bias which arises because the bulk 
flow amplitude is the square root of the quadrature sum of vector 
components and so is a positive definite quantity. Even if the true 
bulk flow were zero, random errors in the components would yield 
a positive amplitude bulk flow. For the same reasons, in the more 
general case of a non-zero bulk flow, random errors bias the ampli- 
tude high. The "error-bias corrected" value of the SMAC bulk flow 
is 620 km s _ . 

Fig.Qshows a "tadpole diagram" of the SMAC data projected 
onto a plane in which the negative X-axis is along the direction 
of the SMAC bulk flow and the vertical axis points to the Galac- 
tic poles. Note the excess of outflowing clusters on the left-hand 
side and the inflowing clusters on the right-hand side, the signature 
of a bulk flow. In contrast, the objects at the top and bottom scat- 
ter around zero peculiar velocity. The bulk flow has a negligible 
component in the vertical (NGP/SGP) direction. Fig. [2] shows the 
SMAC data in four different planes rotated by 45°around the SGY 
axis. Panel (a) corresponds to the Supergalactic Plane. Panel (b) is 
closest to the plane of the SMAC flow. 

It is difficult to illustrate the peculiar velocity errors in tadpole 
diagrams. Fig. [3] shows four Monte Carlo realizations of the pecu- 
liar velocity field where we perturb the observed distances around 
their measured values with a Gaussian random error. Robust re- 
gions of the peculiar velocity field include the clump of positive 
peculiar velocities in the top left quadrant and the negative peculiar 
velocities in the distant part of the bottom right quadrant. 

The bulk flow errors quoted above are marginal errors, i.e. they 
are the diagonal elements of the covariance matrix. The error co- 
variance matrix is not isotropic; it is a triaxial ellipsoid. The long 
axis of the error ellipsoid is oriented toward I — 223°, b = —10° 
(and its antipode). The error along this direction is 222 km s _1 . The 
intermediate axis is toward I = 332°, b = —41°, and the error is 
152 km s _1 . The error along the short axis is 102 km s _1 . The cor- 
responding direction (toward / = 312°, b — 46°) is the direction 
along which the bulk flow is most precisely measured. Because the 
SMAC sample has good sky coverage, the bulk flow is almost inde- 
pendent of the monopole (AH/ H) term: the correlation coefficient 
between it and the bulk flow amplitude is only 0.09. 



We find that the recovered SMAC bulk flow depends very 
weakly on the value of a t h- If we fit for <r t h we find a t h = 
190 ± 105 km s" 1 . For a th = 190 km s"\ the bulk flow is only 
16 km s _1 lower than for the default a t h = 250 km s _ case. 

The flow therefore appears to be quite cold, i.e. the bulk flow 
dominates over sma ll-scale "thermal" moti ons of clusters. The cos- 
mic Mach number l Ostrike r & Sutd ll9901 is the ratio of the bulk 
flow velocity of a volume to the 3D velocity dispersion of the ob- 
jects in the frame of the mean motion, M = j V \ / ( v3 <7 t h ) ■ For the 
SMAC sample, we find M = 2.0 ±1.3 The errors on this quantity 
are obtained by propagation of errors, a procedure not strictly valid 
given the size of the errors in comparison to the measurement. 

To test the coherence of the bulk flow, we have divided the 
sample into statistically-independent nearby (z c < 7500 km s" 1 , 
with an effective depth i?» = 50/i~ 1 Mpc) and a distant subsamples 
(z c > 7500 km s~\ with R t = 98 /i" 1 Mpc). The nearby sample 
has a bulk flow of 640 ±290 kms" 1 toward / = 266°, b = -17°. 
This is significantly different from at the 97% confidence level. At 
thi s depth, flows of order several hundred km s _1 are known to ex- 
ist ICourteau & Dekel200lUZaroubilE o02^. and the SMAC result 
is consistent with these. The distant subsample yields a somewhat 
larger bulk flow (975±300km s" 1 toward I = 267°, b = 16°). Al- 
though the errors are large, the bulk flow of the distant subsample is 
significantly different from zero at the 99.5% confidence level It is 
not statistically different from the bulk flow of the nearby subsam- 
ple. Fig. [4] shows the components of the bulk flow for nearby and 
distant subsamples separated at different values of R. At no value 
of R are the bulk flows of the two subsamples inconsistent. This co- 
herence suggests that there are significant contributions to peculiar 
velocities in the nearby Universe arising from density fluctuations 
on very large (~ 100ft -1 Mpc) scales. 



3.3 Possible Sources of Systematic Error 

In this section, we examine possible systematics and assess the ro- 
bustness of the bulk flow. There are many sources of error that could 
affect a cluster's peculiar velocity, but in general such errors will 
only increase the scatter in the FP or add a random error to the 
cluster peculiar velocities. The bulk flow is the dipole of the pecu- 
liar velocity field, so in order for it to be affected by a systematic 
error, the systematic error must be coherent over large areas of sky. 
The most likely sources of such an effect are problems matching 
data from different observing runs, or systematic errors in Galactic 
extinction corrections. 



3.3.1 Spectroscopic Data 

Systematic errors in velocity dispersion measurements could, in 
principle, lead to large errors in the bulk flow. For example, if ve- 
locity dispersions measured in the North were systematically 1% 
smaller than those measured in the South, then distances in the 
North would be too short by 1.4%. For the SMAC sample, this 
would yield a spurious bulk flow of ~ 50 km s _1 toward the North 
Celestial Pole. Because of these possible effects, great care was 
taken to obtain comparison data for the same objects from multiple 
telescopes (SMAC-I). In SMAC-III, we calculated the corrections 
required for each run as well as the correlated uncertainties in these 
corrections. The latter are fully propagated to the bulk flow and 
are thus included in the error covariance matrix appearing in Equa- 
tion ^2j. Because of the large number of comparison data, these 
errors do not dominate the error in the bulk flow: had we neglected 
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Figure 1. Peculiar velocities of SMAC clusters projected onto a plane in which the negative X-axis is along the direction of the SMAC bulk flow and the 
vertical axis points to the Galactic poles. Clusters within ±45° of the plane are plotted. The circle indicates the estimated distance to the cluster, projected 
onto the plane, and the end of the tail is at the CMB-frame cz, so the peculiar velocity is indicated in km s _1 by the length of the tail. Clusters with peculiar 
velocities away from the origin are filled circles with solid tails, whereas clusters with peculiar velocities toward the origin are shown as open circles and 
dotted tails. The size of the circle scales inversely with the distance errors. The hatched region indicates the Zone of Avoidance < 20°). 



these, we would have quoted an error of 188 km s _1 (compare with 
the correct value of 203 km s _1 ). Subtracting these in quadrature, 
we can quantify the error in the bulk flow from system-matching 
uncertainties: 76 km s~ . 

In order to examine in more detail the effects of individual 
spectroscopic datasets, we have performed a "jackknife" test in 
which we exclude each dataset in turn, and recalculate the bulk 
flow. Some clusters, however, were observed only in one run, so 
the jackknife removes entire clusters from the sample and conse- 



quently changes the spatial sampling of the survey. The results of 
this test are shown in Fig.|5| No single spectroscopic dataset has an 
effect on the bulk flow at a level of more than 70 km s _1 . 



We conclude that errors in spectroscopic systems are fully 
quantified and controlled in the SMAC sample. 



Streaming motions of galaxy clusters — V. The peculiar velocity field 5 




Figure 2. Peculiar velocities of SMAC clusters with symbols as in Fig.0 Four planes are shown: (a) the Supergalactic Plane; (b) a plane rotated by 45° from 
the Supergalactic plane around the the SGY axis. The horizontal axis is 1/a/2 SGX +1/V2 SGZ (toward SGL= 0,SGB= 45°). The SMAC flow is only 
13° from this axis. The signature of a bulk flow, namely infalling objects on one half of the plot and outflowing ones on the other, is clearly seen in this panel, 
c) The SGX=0 plane d) a plane rotated a further 45° so that the horizontal axis is l/\/2 SGX — 1/\/2 SGZ. Only clusters within ±22.5° of the plane are 
plotted. 



3.3.2 Extinction Corrections 

The proximity of the SMAC flow to the Galactic plane suggests 
that errors in the extinction corrections could affect the bulk flow. 
In this subsection, we investigate the effects of extinction on the 
SMAC sample. Note that errors in distance due to random errors 
in extinction are included in the peculiar velocity errors tabulated 
in SMAC-IV, by assuming that the SFD extinctions are accurate to 
16%. This error is applied to the cluster as whole rather than to indi- 



vidual galaxies within the cluster, but is assumed to be independent 
from cluster to cluster. Extinction errors would have the greatest 
effect on clusters with the highest extinction. In the SMAC sample, 
there are 7 clusters (A0400, A0539, A0426, J8, A2657, A3526, 
A3733) with mean E(B - V) > 0.1 mag. and hence A R > 0.26 
mag. We find, however, that when these 7 clusters are excluded, the 
bulk flow drops by only 35 km s _1 . Thus the SMAC bulk flow is 
not substantially affected by random extinction errors. 

Of greater concern is the possibility of coherent errors in the 
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Figure 3. As in Fig. Q] but showing Monte Carlo realizations of the SMAC peculiar velocity field where we have perturbed the observed distances around 
their measured values according to their errors. Notice that the outflow pattern on the left hand side and the inflow pattern in the lower right corner are robust 
to random errors. 



SFD maps. iHudsori 1 19991) tested the SFD maps using data for 
early-type ga laxies, which have a tight intrinsic (B — V) - Mg2 
relation, from F aber et alJ i!989). He found no strong evidence for 
systematic dipolar errors in the SFD maps, and set an upper limit 
of 16% to such systematic errors. lBlakesleeet"aT]<200ll) . using the 
( V— /)-Mg2 relation also found no evidence of a systematic dipole 
error. 

As a further test of systematic extinction errors, we have 
re-computed peculiar velocities and bulk flows using the extinc- 
tion maps of iBurstein & Heiiesl 1:1982, hereafter BH) in place of 
those of SFD. It is worth noting that the BH corrections yield 



some anomalously large peculiar velocities. For example, A2634 
has a peculiar velocity of —1200 ± 390 kms -1 with the BH 
corrections, whereas with the SFD corrections its peculiar veloc- 
ity is —710 ± 370 kms -1 . Similarly, the peculiar velocity of 
A2657 is -3200 ± 1200 kms -1 with the BH corrections, and 
-1600 ± 1200 km s" 1 with the SFD corrections. Using the BH 
corrections, we obtain a bulk flow of 592 km s _1 , a reduction of 
95 km s _1 compared with the result obtained with the SFD correc- 
tions. Thus uncertainties in Galactic extinction appear to affect the 
bulk flow at a level less than half the random error. 
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Figure 4. Bulk flow for different subsamples of the SMAC data selected 
by distance. Solid curves indicate subsamples with data in an inner sphere 
extending from the LG out to R. Dotted curves show bulk flows of subsam- 
ples in an outer shell extending from radius R to the edge of the SMAC 
sample. For clarity, the curves are offset slightly in the horizontal direction. 
For a given curve, the data points at different R are not independent, but at 
each R, the two subsamples are statistically independent. In all cases, the 
bulk flows of inner and outer samples are consistent with each other. 



Figure 5. The jackknife test applied to the bulk flow statistic showing the 
effect of excluding a given spectroscopic run from the full solution. The 
vertical axis shows the amplitude of the bulk flow and the horizontal axis in- 
dexes different spectroscopic runs. Open circles show results for each spec- 
troscopic run in turn. Those runs whose exclusion changes the bulk flow by 
more than 50 km s -1 are labelled. For details of the labels, see lSmith et alj 
(2000). The solid line shows the best fit result from the whole sample, and 
the dashed lines show the 1 a errors on the amplitude. No individual spec- 
troscopic run has a significant effect on the bulk flow. 



3.3.3 Stellar Populations 

Variations in age and metallicity will introduce extra scatter in the 
FP. If there are systematic differences from cluster-to-cluster, then 
it is possible that the cluster peculiar velocities could be signifi- 
cantly in error. Mg2 can serve as an indica tor of age and/or metal- 
licity. For example, the models of Worthev 1 1994) indicate that, for 
a typical elliptical, a —0.15 change in metallicity at solar abun- 
dance yields a Mg2 change of —0.03 and an R-band change of 
—0.13 mag. Similarly if we add a 0. 10 mass fraction in the form of 
an intermediate-age (5 Gyr) population, Mg2 changes by —0.0083 
mag. and R-band light changes by —0.10 mag. 

In SMAC-IV, we compared residuals from the Mg2-<r relation 
with residuals from the inverse FP on a cluster-by-cluster basis and 
found no strong evidence for a correlation. 

We also showed, however, that the FP-Mg2 relation, which in- 
cludes Mg2 as a parameter in addition to the usual 3 FP parameters, 
does reduce the scatter in log a in the inverse FP. We prefer not to 
use the FP-Mg2 relation for our default solution for two reasons: 
first, Mg2 is not available for all of our galaxies, and second, the 
scatter in distance for the FP-Mg2 relation is actually larger than 
for the inverse FP, because the log R c coefficient in the FP-Mg2 
relation differs from that of the FP relation. Nevertheless, if there 
are systematic variations from cluster-to-cluster, these can be cor- 
rected by including Mg2 in the distance indicator. However, when 
we compare the same sample of clusters, we find that including the 
Mg2 term increases the bulk flow by 261 km s _1 to 831 km s~ 
towards lb254-7, while also increasing the error in the bulk flow to 
230km s _1 . 



3.3.4 Morphological Mix 

The SMAC sampl e contains both E and SO galaxies. In 
iHudson et alJ 1 19971 hereafter PP-II), we compared the FP rela- 
tions of E (T -4) to SO (T > -3) types. There we found a 
small and marginally significant (2<j) offset in the zero-point, in the 
sense that E types are observed to hav e larger velocity dispe rsions 
than SO galaxies at fixed R c and (fi) c . Blakeslee et al. (2002J) also 
noted differences between E and SO in their comparison of FP and 
Surface Brightness Fluctuation (SBF) distances. In SMAC-IV, we 
repeated this analysis, this time with a free coefficient of the RC3 
T-type. There we found a small and marginally significant correc- 
tion to the FP predicted velocity dispersions (—0.0043 ± 0.0019 
in log a per unit of T-type). This corresponds to E types having 
2.7 ± 1.3% larger velocity dispersions at a given R c and (^i) o than 
SO galaxies. This difference is smaller than, but consistent with, the 
offset found in PP-II. We have not included this term in our default 
solution because we do not expect any significant difference in the 
ratio of Es to SOs across the sky. When the term is included, the bulk 
flow drops by only ~ 30 km s _ . If we treat E and SO subsamples 
separately, they both give a consistent bulk flow, the E sample of 
417 galaxies yields a slightly lower amplitude (601 ±195 km s _1 ) 
bulk flow than the SO sample of 277 galaxies (762 ± 294 km s^ 1 ). 
The larger bulk flow error for the SO-only sample arises partly from 
the smaller number of galaxies in the SO sample, and partly from 
the larger scatter in the SO FP relation. 
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Figure 6. As in Fig.|5] but showing the effect of excluding individual clus- 
ters of galaxies. 



FP distances (as opposed to its redshift), and so these distances 
must be corrected for Malmquist bias. One advantage of working 
with clusters is that this correction is small, since it scales as the 
inverse of the number of objects in the cluster. 

It is also possible to fit the flow model and FP relation simul- 
taneously. In this alternative "Method II" fit, redshift (corrected by 
the flow model) is the a priori distance indicator. As a result, no 
Malmquist bias are corrections necessary 1 . A Method II analysis 
of the SMAC sample assuming a simple bulk flow model yields a 
bulk flow amplitude only ~ 20 km s~ less than our standard re- 
sult with Malmquist bias corrections, indicating that the details of 
Malmquist bias correction have little effect on the bulk flow. 

3.3.7 Summary of Systematic Effects 

We have examined several possible sources of systematic error. The 
most obvious source of systematic error, namely mismatches be- 
tween velocity dispersions obtained on different runs is fully ac- 
counted for in our standard error analysis. The only systematic ef- 
fect which might reduce the bulk flow is using the BH extinction 
maps (at the price of introducing very large peculiar velocities for 
two clusters). This would lower the flow by less than 100 km s _1 . 
Thus systematic errors in the bulk flow appear to be negligible in 
comparison to the random error of ~ 200 km s _1 . 



3.3.5 Outlying Data 

In SMAC-IV, we removed 5 galaxies which deviated by more than 
3 sigma from the best-fitting FP. Re-including these 5 outliers re- 
duces the bulk flow by 56 km s _1 . 

In Fig. |^| we show the effect on the bulk flow amplitude of 
removing entire clusters. Clusters which affect the bulk flow more 
than most are indicated. No single cluster influences the bulk flow 
upwards by more than 55 km/s. In particular, removing the nearby 
Centaurus cluster (which includes both the Cen30 and Cen45 sub- 
components) has little impact on the bulk flow (reducing it by only 
40k mr'). 

Gibbo ns"et alj 1200 lh have suggested that FP clusters with 
large internal FP scatter have significantly larger peculiar veloc- 
ity amplitudes. This might arise, for example, in cases where it is 
difficult to separate double clusters in which the two components 
are superimposed along the line of sight but are at slightly different 
distances. The spread of distances would then be increased, and the 
mean distance, and hence peculiar velocity, would be more suscep- 
tible to how the two components are sampled. Classic cases where 
this situation occurs are Abell 400 and Cen30/45. We examined 
this issue in Paper IV and found no strong evidence for a difference 
between the peculiar velocities of high and low scatter clusters. 

Despite our non-detection of this effect, we have experimented 
with a weighting scheme in which the error for each cluster's dis- 
tance is based on the FP scatter for that cluster, rather than a global 
value. This procedure down-weights those clusters with large inter- 
nal scatter. When we apply this weighting scheme, we find that the 
bulk flow drops, but only by ~ 80 km s . 

3.3.6 Malmquist Bias Corrections 

The procedure we have followed involves first calculating cluster 
distances from the inverse FP and then fitting these to a flow model. 
In the terminology of Strauss & WillicM Jl 995h . this is "Method I". 
In this method, the estimate of each cluster's distance is its inverse 



4 OTHER FLOW MODELS 

4.1 Hubble Bubble 

IZehavi et al] Jl998l) . in an analysis of 44 SNIa distances, claimed 
evidence for a "Hubble bubble", i.e. that the local Universe within 
7000 km s" 1 is underdense with the result that the local value Hub- 
ble constant is to high compared t o the global value, AH/H = 
6.5 ± 2.2%. lGiovanelli et alJil99sl) . using the TF relation, found a 
statistically insignificant "bubble" of AH/H = 1 ± 2%. 

We have tested this model with the SMAC data. We measure a 
local Hubble bubble of 2.3 ± 1.9% within the same distance stud- 
ied by Zehavi et al. This result is consistent with zero. However, 
allowing for errors in the Zehavi et al. result as well as ours, the 
two results are consistent. We do not find any evidence of a Hub- 
ble bubble if we cut the SMAC sample at other radii. Taking all of 
these results at face value suggests there is at most a modest Hubble 
bubble with AH/H ~ 3 ± 1.3%. 

4.2 Bulk Flow plus Shear 

The large scale shear of the peculiar velocity field offers an oppor- 
tunity to identify the distances of the sources responsible for the 
bulk flow iLilie et al. 1986) particularly if they are near the sur- 
vey limits or in the Zone of Avoidance. For example, if there were 
an attractor just beyond the SMAC survey limits, then, relative to a 
pure bulk flow, the peculiar velocities would be more positive at the 
edge close to the attractor and less negative on the opposite side of 
the sky. Thus the residuals from the constant bulk flow would have 
a tidal pattern, which might be measureable. Since the strength of 
the tidal field is inversely proportional to distance cubed, then if the 
attractor was very far distant, the measured shear would be small. 

1 Technically, there is a Malmquist-like bias associated with the error in 
the estimated redshift, but at the distances considered in this paper, this bias 
is very small. 
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To measure the shear, we modify equation Q as follows: 

U(r) = V r + r h r (4) 

where h is a symmetric shear tensor. The shear tensor measures 
anisotropics in the Hubble expansion. When the shear tensor is di- 
agonalized, each eigenvalue corresponds to the difference in the 
expansion rate (AH/H) along the corresponding eigenaxis. This 
differs from equation Q where the Hubble constant was allowed 
to float, but was forced to be the same in all directions. 

When this fit is performed, we find no significant reduction in 
\ 2 compared with bulk flow only fits. The direction of expansion 
(6 ± 3%) is towards I = 315°, b = -7° and its antipole. The di- 
rection of expansion is very close to the negative SGY axis in the 
Supergalactic Plane. The direction of compression (5.5 ± 3%) is 
towards I = 208°, b — —68° and antipole. However, these ampli- 
tudes are not statistically significant, and consequently the direc- 
tions of expansion and compression are ill-defined due to the large 
errors. With this fit, the bulk flow is essentially unchanged: we find 
615 ± 211km s" 1 toward I = 259°, b = 4°. 



4.3 Attractor Models 

Next we consider models with simple attractors. Although such 
models are likely to be only a crude approximation to a Gaussian 
random density field, they allow us to gain some insight into cos- 
mographical features in the nearby Universe. The SMAC cluster 
data are rather sparse so we do not attempt to identify attractors 
in the data themselves. Instead we concentrate on attractor models 
using distances and profiles fixed by other authors. 

We model the radial infall component towards an attractor 
with the following functional form: 



Ma(r) 
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(5) 



where d a is the distance from the LG to the attractor, r a = d a — r 
is the vector from the point r to the centre of the attractor, c a is 
a core radius and v n is the velocity infall towards the attractor at 
the position of the LG. The infall is it a ~ at r a ~ 0, peaks near 
r a ~ c a and then falls of as u a ~ r^"™" at large r a . Here we fix 
n a = 2. 

The first attractor we consider is the "Great Attractor", us- 
ing the flow model of Faber and Burstein (1988, hereafter FB88), 
except that for simplicity we use n a = 2 rather than n a = 1.7 
adopted by FB88. We keep the location (t-qa = 4200 km s" 1 in 
the direction I = 309°, b = 18°) and core radius (1400 km s _1 ) of 
this attractor fixed with the parameters in the FB88 model, and fit 
only for one free parameter: the infall towards the GA at distance 
of the LG. With no bulk flow in the fit, the infall from the GA is 
88 ± 68 km s _1 . This is consistent with zero and differs signif- 
icantly from the result of FB88 who found dga = 535 km s _1 . 
Note that in Fig.|2ji, there is no evidence of infall into the GA on 
its the far side. We caution, however, that our data are sparse in the 
immediate vicinity of the GA. There are only three SMAC clusters 
for which the FB88 GA model predicts significant negative peculiar 
velocities (mga < —500 km s -1 ). For each of these three clusters, 
the SMAC peculiar velocity is positive. If we include a bulk flow in 
the fits, the GA infall is even smaller (13 ± 70 km s -1 ) . These con- 
clusions do not change if we adopt the GA position of Kol att et ail 
( 1995). We conclude that there is little evidence from our sample 
for a very large mass at the position of the GA. However, it should 



be noted that the errors are large: the 95% upper limit from our bulk 
flow plus GA mode l is vga ~ 150 km s _1 . 

Hudson 1 1993), in his analysis of the density field of optically- 
selected galaxies, found the GA overdensity was better described as 
a broad overdensity centred on the Centaurus cluster. He predicted 
an infall of 287 ±62km s" 1 at the L G for /3 opt = 0.5. In their anal- 
ysis of SBF data, Tonrv et al. 1 2000) fitted a "GA" model, the centre 
of which is almost coincident with Centaurus cluster. They obtained 
an infall of 2 89=1=137 km s" 1 at the LG, in good agreement with the 
Hudsqi] 1 19931) prediction. I f we adopt the GA/Centaurus attractor 
centre o flTonrv et alj iSoOO). we find an infall of 58 ± 153 km s _1 . 
It is difficult to ma ke an exact comparison because the details of 
the infall model of lTonrv et alJ {2000) is different from ours, but 
our results do not appear to be in conflict with theirs. 

It is clear that local (R < 60/t~ 1 Mpc) attractors are not 
responsible for the large-scale flow seen in the SMAC sample. 
The sample must be responding to sources at distances beyond the 
SMAC effective depth of 6000 km s~ . The next attractor consid- 
ered is the Shapley Concentration (SC) centered on the rich cluster 
Abell 3558 at I = 312°, b = 31° and a distance of 145 /i _1 Mpc. 
Since there are SMAC clusters quite close to the SC, it is important 
to model the SC as an extended mass distribution. The model given 
by equation ^5) allows for a core radius of the attractor. The SC 
core radius is not well determined by the SMAC data themselves (a 
fit yields ~ 30 ± 20h^ Mpc) so we fix it at 30 /i _1 Mpc. 

If we fit for only a SC attractor with no bulk flow, we find an 
infall at the LG of 200 ± 60 km s _1 . This corresponds to an excess 
mass of 9 ± 3x 10 16 (S] m /O.3)°' 4 /i" 1 M from the SC region. 

If we fit for the SC attractor plus a bulk flow, the infall due 
to the SC reduces to 140 ± 80 km s _1 , and so is not statistically 
significant. For this model, the bulk flow is 620 ± 220 km s _1 and 
the direction has swung around by ~ 10 deg to I = 248°, b — —6°. 
This bulk flow is inconsistent with zero at the 95% CL. Note that 
there is strong covariance between the SC infall and the bulk flow 
which boosts both of their marginalized errors. 

We conclude that there is a tantalizing suggestion for substan- 
tial mass at the SC, but that it is not the sole source of the SMAC 
bulk velocity. 



5 THE PECULIAR VELOCITY FIELD PREDICTED BY 
THE IRAS PSCZ SURVEY 

5.1 Method 

The toy attractor models discussed above are rather crude. A better 
approach to modeling the peculiar velocity field is to use a red- 
shift survey to reconstruct the real space density field of galaxies. 
If mass density contrasts are related to galaxy number-density con- 
trasts according to a simple biasing relation, e.g. 5 e = bS, then we 
can obtain predictions for the peculiar velocity field using the linear 
theory equation 



JL 

47T 



r' 



(6) 



where (3 = f2 m °' 6 /6. This method has been used by a number 
of workers to compare observed and predicted peculiar veloci- 
ties and hence obtain (3 ( see r eviews by j^jraus s & Willickl ll995h . 
ICourteau & Dekeil i2001) and Zaroubif |2002|) ). For predictions of 
the peculiar velocities of clusters of galaxies, linear theory is ade- 
quate. 

While the integral in equation {6) extends over all space, in 
practice, redshift surveys typically do not have data in the Zone of 
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Avoidance and are truncated at large distances where the correc- 
tions for selection effects become large. In the case of a distance- 
limited redshift survey, equation then yields peculiar velocities 
in the frame of the center of mass of the redshift survey volume. 
In general, there will be contributions to the CMB -frame peculiar 
velocity from sources outside the redshift survey volume. If the 
redshift survey is much deeper than the peculiar velocity survey, 
and there are no massive structures in the ZOA, then the tidal or 
quadrupole effect from distant attractors, which falls as r~ s , will 
be small. It is then sufficient to model these residual contributions 
by adding a simple dipole term V cx t to equation JSJ: 



u(r) = ui hl (r) + V c 



(7) 



There are several approaches to fitting (3. One common ap- 
proach is to perform the comparison between predicted and ob- 
served peculiar velocities in the LG frame. In this comparison, 
the exterior dipole cancels from the predictions since ulg(V) = 
u(r) — u(0). An alternative approach is to fit in the CMB frame, 
but this ignores the LG as a data point entirely. 

In this paper, we make the comparison in the CMB frame but 
we force our solutions to be consistent with the LG's peculiar ve- 
locity. We implement this by including the three Cartesian compo- 
nents of the LG's peculiar velocity as if they were three additional 
peculiar velocity data with zero observational error. We do, how- 
ever, allow both the LG peculiar velocity components and the pe- 
culiar velocities of the SMAC clusters to have a random "thermal" 
component to their error. In this context, the "thermal" component 
represents the combined effects of non-linearities in the peculiar 
velocity field and errors in the peculiar velocities predicted from 
the galaxy density field. 
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Figure 7. A comparison of observed peculiar velocities and predictions 
from the PSCz density field with ft = 1. SMAC clusters are indicated by 
filled circles. Clusters which deviate by more than 2a from the model pre- 
dictions have error bars plotted and are labelled. Vectorial components of 
the LG's peculiar velocity are indicated by open squares and are labelled. 
The best fitting external bulk flow V cx t has been subtracted from the ob- 
servations. Symbol sizes scales inversely with observational error. The solid 
line shows the ft of best fit; dashed lines indicate the errors. 



5.2 Application to PSCz 

For the galaxy redshift survey, we have used the IRAS Point Source 
Catalogue Redshift (PSCz) surve y fcaunders et all2000bh with pe- 
culiar velocity predictions from Branchini et al. 1 1999). The dis- 
tance limit of the density field is 20000 km/s. The Galactic Plane 
(|6| < 8°) is filled by interpolating the data at higher |6|. The "real 
space" distances are obtained in a self-consistent way by an itera- 
tive method (see lBranchini et all J 19991) for additional details). 

The flow model uses equation Q where Uy m are the IRAS 
predictions, obtained by using the IRAS galaxy density field as 8 g 
in equation JSJ. This model has four free parameters: ft and the 
three components of the V C xt- We set the a t h = 150 km s _1 for 
this comparison. 

One final correction is necessary. We are making predictions 
for clusters of galaxies, which are significant mass concentrations. 
In the IRAS density field, these clusters will be located at a given 
"real-space" position as obtained by the self-consistent iterative re- 
construction method. Due to distance errors, the observed SMAC 
cluster will, in general, be located a different distance along the 
line of sight from the IRAS cluster. We do not want the IRAS clus- 
ter to generate a spurious peculiar velocity at the location of the 
SMAC cluster, so for each SMAC cluster we identify and temporar- 
ily delete its counterpart in the IRAS density field before calculating 
the predicted peculiar velocity. 

5.3 Results 

Simultaneous fits of the parameters in equation Q Me given in 
TableQ The best fit has ft = 0.39 ± 0.17 and V cxt = 372 ± 



127 km s" 1 towards I = 273° ± 17°, b = 6° ± 15°. Fig.0shows 
uun (with ft = 1) vs. iiobs- The slope of the line of best fit then 
gives ft. 

For our best fit, \ 2 = 57.5 for 54 degrees of freedom. In 
comparison, the null case of ft = 0, Vcxt = has x = 95.1 for 
58 degrees of freedom. The marginalized errors quoted above are 
large because there is strong correlation between the ft and V cxt 
terms: as ft increases, |V e xt| drops. The error contours associated 
with this fit are shown in Fig. [8] The covariance complicates the 
interpretation of these results. Of the two components of the flow 
model, ft is the more significant. If V cx t is set to zero, ft = 0.70± 
0.13, a > 5a detection of ft This "high-/?" model is disfavoured 
at the 98% CL compared to the model with ft and V cx t free, for 
which x 2 falls by 10 with the removal of 3 degrees of freedom. 
Finally, if ft is fixed at 0.5 (as external comparisons suggest, see 
below), then we find V ox t = 327 ± 104 km s" 1 . 

It is well known that IRAS misses early-type galaxies in the 
cores of clusters. This effect may be particularly dramatic in regions 
such as the Shapley Concentration. We have fit a flow model with 
an attractor at the position of Shapley (core radius 3000 km s _1 ) 
in addition to the IRAS predictions and an external bulk flow V ex t ■ 
For this fit, we find an infall of 129 ± 84 km s -1 at the LG. This is 
only ~ 10 km s _1 less than the infall into SC found above (without 
the PSCz flow model), indicating that SC is very weak in the PSCz. 
For this fit the external bulk flow reduces to 300 ± 140 km s - . 
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Table 1. Fits using the PSCz predicted peculiar velocity field 
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Figure 8. Parameter fits to the SMAC+LG peculiar velocity data using the 
PSCz predictions. Contours show the 68% and 95% confidence contours on 
01 and the {/-component of V cx t ■ The open circle shows the best fit values. 
The hashed vertical bar indicates the concordance value of /3t discussed in 
Section l5TT1 



(ii) Fitting peculiar velocity data. IWillick & Strauss! < 19981) 
found 0i = 0.50 ± 0.04 by comparing the predic tions of th e 
IRAS 1.2 Jy survey t o the Mark III data JWillick et all 1 1997b . 
iBranchini et id] ll999l) found 0i = 0.4 2 ± 0.04 from the com - 
parison of the PSCz t o the SFI catalogue Giov anelli et alj ll998t) . 
Blakeslee et al. 1 1999) obtained 0i = 0.42±g;J§ from a compari- 
son of Surface Brightness Fluctuation survey peculiar velocities to 
IRAS 1.2 Jy predictions. These three analy ses are conducted in the 
LG frame and are insensitive to V^ t . Izaroubi et a l] |2002|) found 
0i = 0.51 ± 0.06 from comparison of ENEAR Ida Costa et al] 
2000b) and SFI peculiar velocities to PSCz in the CMB fram e. 

(iii) Redshift-space distortions. Hamilton et alj |2000) found 
01 = 0.41 

l"o'.i2 from an analysis o f redshift-space distortions in 
the PSCz survey. iTavlor et alJfcOOll) found essentially identical re- 
sults. 

(iv) Expectations from cosmological models. Usin g the 
WMAP, 2dF and Lyman-a forest data lspergel et alj J20031) found 
O m °' 6 crg = 0.38^0 05- Combining this with the m easurement 
ff 8 ,i = 0.80 ± 0.05 from lHamilton & Tegmarkl <2002h . we obtain 
01 = 0.48 ± 0.06. 

The "concordance" value 0\ — 0.5 is consistent with all of the 
above determinations at the 2a level. Note that the value 0i = 0.5 
requires a residual motion of ~ 160 km s _1 towards / = 319°, 
6 = 13° (or V x = 119, V y = -110, V z = 35) from outside the 
PSCz volume (i.e. from beyond 200 ft~ 1 Mpc or in the Zone of 
Avoidance) in order to match the observed peculiar velocity of the 
LG. 



5.4 Discussion 

The degeneracy between and V cx t make it difficult to build a 
unique flow model from SMAC+LG data alone. There are several 
independent methods of determining 0i and V ex t- 



5.4.1 External constraints on 0\ 

(i) Matching the dipole to the LG's motion. This method con- 
sists of using equation JSJ to predict v(0) and adjusting 0i to obtain 
the best fit with the observed 627 km s _1 motion of the LG. If the 
flow model allows for V e *t as well, then the two components are 
degenerate (one is fitting 3 data with 4 free parameters). As a result 
for this method to yield constraints on 0x one has to make some 
assumptions for V ex t- 

One approach is to allow cosmological models to predict the co- 
herence between a distance-limited gravity dipole and the peculiar 
velocit y of the LG. The most rec ent estimate of 0i by this method 
JCiecielag & Chodorowskil2004l) yields 0i = 0.64+°; 2 i. 



5.4.2 External constraints on V cx t 

What is the expected value of V cx t? There are four sources of un- 
certainty in the PSCz predicted peculiar velocities. 

(i) Shot noise with in 7? ~ 200/i _1 Mpc. At R = 150/i _1 Mpc, 
ISchmoldt etal] 1 1999) quote errors of 160/3i km s 1 on the pre- 
dicted velocity of the LG. They do not extend their analysis to 
greater distances, but an eyeball extrapolation of their Figure 5 sug- 
gests that this grows to ~ 200/3i at R = 200fe~ 1 Mpc 

(ii) Sources in the Galactic Plane. Saunders et al. 1 2000a) pub- 
lished preliminary results from the "Behind the Plane" (BTP) sur- 
vey (an extension of the PSCz to lower Galactic latitudes). The 
BTP dipole in units of km s~ (assuming 0i — 0.5), grows from 
V x = -60, V y = -325, V z = 350 at R ~ lOO/iT 1 Mpc, to V x = 
-60, V y = -475, V z = 350 at R ~ 200/T 1 Mpc. The misalign- 
ment with the CMB dipole of only ~ 13°. At R ~ SOOh" 1 Mpc, 
the dipole is V x = -25, V y = -675, V z = 350, or ~ 761 km s" 1 
towards I = 268°, b = 27°, and the misalignment is reduced to 
~ 8°. The BTP dipole is in better directional agreement with the 
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LG's observed motion than the PSCz dipole. It is possible that this 
alignment is fortuitous since the shot-noise effects (and systemat- 
ics such as the "rocket-effect" JKaiserll987h ') are likely to be large 
beyond 200 ft _1 Mpc. There is a particularly strong change in V y 
in the range 180-240 h~ Mpc which is present in both the PSCz 
and BTP dipoles, but is much stronger in the latter. This suggest 
that there are dynamically important structures close to the Galac- 
tic Plane. 

For purposes of assessing V cx t result from the SMAC sample, 
an estimate of the effect of hidden sources in the Plane is given by 
the difference between the growth in the BTP dipole and the PSCz 
dipole at distances beyond the SMAC sample (~ lOO/i" 1 Mpc). 
The difference in km s" 1 for (3i = 0.5 is A 14 = +105, AV y = 
-135, AV Z = -20, or 172 km s _1 towards I = 308°, b = -7°. 
We estimatethe error on this to be 50% or 85 km s _1 . 

The "Clusters Behind the Zo ne of Avoidance" (CI ZA), an X-ray 
selected survey galaxy clusters, Ebeling et al. 1 2002) also indicates 
several massive clusters at low galactic latitudes. However, these 
appear to be located either in the vicinity of the GA at 60 h~ 1 Mpc 
(A3627 and CIZAJ1324.7-5736) or nearer to SC at a distance of 
~ 150/i -1 Mpc (CIZA J1653.0-5943 and Triangulum Australis) 
and not at 180-240 /i _1 Mpc suggest ed by the BTP dipol e. The 
gravity dipole of the CIZA sample iKocevski et alj Eo03) has a 
strong contribution at ISOh' 1 Mpc which is not seen in the PSCz 
or the BTP. 

(iii) Sources beyond 200 ft, _1 Mpc. Since the PSCz/BTP data 
are noisy beyond 200 h~ 1 Mpc, it is difficult to assess this empir- 
ically. For a ACDM cosmological model with T = 0.21, fl m = 
0.3, n = 1 and erg = 1, we expect the rms contribution to the bulk 
motion in the nearby Universe arising from sources beyond 200 
h~ 1 Mpc to be ~ 50 km s _1 in each component. 

(iv) Missing contributions from high-density superclusters. 
IRAS does not detect early-type galaxies. For individual clusters 
this is a weak effect, which can be compensated by increasing 
(3\. There are certain extreme regions, such as the SC and the 
Horologium-Reticulum superclusters, where there is an astound- 
ing excess of clusters. For example, iTullv et alJ <1992h find 29 
Abell/ACO clusters within 50 ri -1 Mpc of the centre of the SC. This 
corresponds to a mean overdensity S c ~ 5 in clusters. If 6 c /6i ~ 4, 
(Branchini et al. 1999), then we expect the PSCz overdensity to be 
Si ~ 1.25, yielding an infall at the LG of 250/3i km s _1 . In fact, 
the observed PSCz overdensity is only 0.2, yielding an infall of 
only 40/3i km s -1 . A similar situation occurs for the Horologium- 
Reticulum supercluster. Tully et al. finds a cluster overdensity of 
~ 4 on scales of 50 /i _1 Mpc but the PSCz overdensity is only 0. 14. 
Clearly linear biasing does not operate in these very high-density 
regimes. 

Because the cluster overdensity is likely an overestimate of the 
true mass density (even when scaled by bi/b c ~ 0.25), and the 
PSCz overdensity an underestimate, one might expect the true sit- 
uation to lie in between. For /3i = 0.5, we will estimate a residual 
dipole of ~ 50±25 km s _1 in the direction of the SC, and, because 
it is at greater distance, half of that for Horologium-Reticulum. 

Adding the BTP residual dipole (172 km s _1 toward 
I = 308°, b = b = -7°), plus contributions directed to SC (50 
km s" 1 ) and HR (25 km s" 1 ), we obtain an estimate of \^ e xt — 
225 km s -1 towards I = 305°, b = -4°. Most of this is due to the 
extra Galactic Plane sources in the BTP dipole. Errors on V cyi t esti- 
mated to be 90 km s" 1 (systematic uncertainty) and 75 (shot noise) 
and 50 (residual beyond 200 h~ 1 Mpc). Summed in quadrature, 
this yields an error of 150 km s _1 . Clearly, given the uncertainties, 



there is no firm external evidence that V cyi t is required. The fact 
that the best estimate is an agreement with the measured direction 
and amplitude of V cxt from the SMAC survey (372 ± 127 km s _1 
towards I = 273° ± 17°, b = 6° ± 15°) suggests that these exter- 
nal sources do exist. 



6 COMPARISON WITH OTHER RESULTS 

The bulk flow calculated in this paper is a statistic derived from the 
SMAC sample. This is not to suggest that the true flow field within 
12000 kms -1 is a pure dipole bulk flow. Indeed, a pure dipole 
velocity field is unlikely given that the mass density field is a Gaus- 
sian random field with fluctuation power on a range of scales. Due 
to sparse sampling, the bulk flow of any peculiar velocity sample 
will have contributions from the true bulk flow of the volume, as 
well as from higher order multipoles. This effect was first empha- 
sized by Watkins & Feldma^ ll995t) . It would, therefore, be naive 
to compare the bulk flow from different sparsely-sampled surveys 
without allowing for the error introduced by sparse sampling. Nev- 
ertheless, that has been the approach taken in the past, and it has 
led to an exaggeration of the differences between various studies. 

A full statistical comparison between surveys, based on the 
three components of the bulk flow and the full error covariance ma- 
trices allowing for the above-mentioned sparse sampling effects, is 
beyond the scope of the present paper. A preliminary comparison 
of the bulk flows of SMAC lWillicM ll999l) . LP. lDale et al1ll999l) . 
Colless et al. lereafter EFAR), and the SNIa compilation of 

Tonr^HrtT 20031) has been performed by iHudsonl 1200 3l) . who 
concluded that there was no inconsistency between any of these 
surveys, except possibly for LP which is inconsistent at the 94% 
level. 

Instead of the full vectorial comparison, in this paper we opt 
for a simpler, illustrative comparison which considers simply the 
component of the bulk flow projected along a given axis directed 
towards / = 300°, b — 10°. This "concordance" direction has been 
chosen because i t is close to the mean flow of the above-mentioned 
sparse surveys lHudsorl2003h . is within a few degrees of the Super- 
galactic Plane, and is also ve ry close to the bulk motion originally 
reported for the 7S sample by Dressier et al. 1 1987). In his analysis 
of the Mark II peculiar velocity catalogue, Hudson 1 1994) found 
that the gravity of the local density field to 8000 km s _1 could not 
account for a residual bulk motion of 400 ± 45 km s~ x toward 
I = 292°, b — 7°. If this model is correct, then deeper surveys such 
as SMAC should have bulk flows close to this residual motion. 

For the SMAC sample, with the bulk flow fixed in the con- 
cordance direction / = 300°, b = 10°, the best fitting amplitude 
is 400 ± 120 km s _1 . The following sections compare our results 
with those from from other surveys, focussing first on those based 
on sparse samples on large scales . 

6.1 Large-scale Surveys 

6.1.1 Brightest Cluster Galaxies 

Lauer & Postma rt Jl994h found a bulk flow of 689 ± 178 km s" 1 
towards / = 343°, b — +52° from 119 clusters using the a dis- 
tance indicator based on the photometry of brightest cluster galax- 
ies (BCGs). The direction of the SMAC bulk flow is nearly 90° 
away from that of LP, so these two results appear to be in poor 
agreement. No r does LP agree w ith results from other peculiar ve- 
locity surveys toaleetalj|l999l EFAR). The LP sample is denser 
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than other all-sky surveys on simuilar scales and so it is less af- 
fected by sampling errors than other sparser surveys, including 
SMAC. Therefore sparse sampling is unlikely to be a source of 
disagreement between LP and other surveys. 

In SMAC-IV, we compared our cluster distances with those 
of LP. We found that of 41 clusters in common, four were dis- 
crepant at greater than the 2a level: A262, A1060 (Hydra), A3381 
and A3733. In all cases, comparison with the observed redshifts 
indicates that the LP distances are too long corresponding to BCG 
magnitudes whch are too faint. 

One possible explanation for th is discrepancy is obscuration 
by dust in BCGs. Laine et al. (2003) examined dust in 75 of LP's 
1 19 BCGs. They found signs of dust in ~ 38% of BCGs. In some 
cases the dust appears to be quite extended, and this could affect 
the LP magnitudes. 

There is evidence that galaxies classified by Laine et al. as 
having filamentary or patchy dust are systematically too faint. If 
we compare the ratios of the estimated distance from the LP BCG 
method to the distance estimated by assuming that the BCG is at 
rest with respect to the CMB, we find that for 49 galaxies with no 
dust the distance ratio is 0.98±0.02. Galaxies with only nuclear (< 
1") dust disks (Laine code "D"), dust spirals (code "S"), or rings 
(code "R") have a distance ratio 0.92 ± 0.04. These small scale 
features are unlikely to affect the 10 /i _1 kpc aperture magnitude 
which correspon ds to a much larger angular scale, typically tens of 
arcseconds (Postman & Lauer 1995). However, the distance ratios 
are higher for galaxies with filamentary ("F") dust (1.06 ± 0.06), 
patchy ("P") dust (1.12 ± 0.08) or both patchy and filamentary 
dust (1.28 ± 0.10). The latter category consists of only 4 galaxies, 
the BCGs of A0262, A1060, A3698 and A3733. Three of these 
(A0262, A1060 and A3733) clearly have very extended dust in gray 
scale images (fig. 1 of Laine et al.). These same three galaxies are 
extreme outliers in the SMAC vs. LP comparison discussed above. 

If we remove from LP's sample the 19 galaxies with dust 
classified as filamentary or patchy (or both) in Laine et al.'s 
Table 2, we obtain a bulk flow of 707 ± 261 km s _1 toward 
I = 336° ± 46°, b = 38° ± 23°. The error-bias-corrected ampli- 
tude is 384 km s _1 (to be compared with LP's quoted 689 km s _1 ). 
A3381 is not in the Laine et al. sample, but is the most dis- 
crepant cluster in the SMAC-LP distance comparison performed 
in Paper IV (3.7 a). While this cluster has a somewhat larger in- 
ternal FP scatter Ijorgensen et aljfl996h . its SMAC peculiar ve- 
locity is not significant 1103 ± 679 km s~ . In contrast, its pe- 
culiar velocity according to LP is large and signficantly differ- 
ent from zero (—5900 ± 2100 kms -1 ), so we suspect that for 
this object LP are in error, rather than SMAC. If we also remove 
A3381, the LP bulk motion becomes 667 ± 242 km s" 1 towards 
I = 315° ± 48°, b = 35° ± 26°. This differs from no bulk flow, 
but only at a marginal (95%) confidence level. This is in better 
agreement (both in direction and amplitude) with the SMAC bulk 
flow. Along the direction I — 300°, b = 10°, the motion of this 
sample is 700 ± 380 km s _1 . Note that because of the error el- 
lipsoid is triaxial it is possible that a flow solution fixed along a 
given direction has a higher amplitude than that of the flow along 
the direction of best fit. 

There is therefore a strong hint that the original LP sample 
was affected by dust. Although it is tempting to use the result from 
the culled sample as a "corrected" LP bulk flow, it is important to 
remember that there remain 25 galaxies in the BCG sample of LP 
which are neither in Laine et al. sample nor have comparison pe- 
culiar velocity data from SMAC. It is likely that some (~ 25%) of 
these will be strongly affected by dust. Furthermore, it is also possi- 



ble that the LP motion is affected by other systematics, such as the 
tendency for brighter BCGs to inhabit more X-ray-luminous clus- 
ters iHudson & Ebelinslll997t) . We conclude that the original LP 
result is in marginal disagreement with SMAC, but this disagree- 
ment is no longer significant once some BCGs with clear evidence 
of dust are removed from the LP sample. 



6.1.2 Tully- Fisher distances to clusters 

iMould et al] 1 19931) observed 38 TF clusters and fitted a bulk flow 
of811±138 kms -1 towards / = 332°, b = -15°. Along the "con- 
cordance" direction, this projects to 617 ± 138 km s~ . Exclud- 
ing two outliers (NGC 5419 and A3627), they obtained a flow of 
559±107kms _1 toward / = 326°, b = -9° or 473±107 km s" 1 
tow ard 1 = 300°, b = 10°. 

IWillickl j 19991) found a bulk flow of 961 ± 280 kms" 1 to- 
ward I = 266°, b = 19° from 15 clusters using the Tully-Fisher 
relation. The amplitude and direction are in excellent agreement 
with SMAC bulk flow, although the er rors are large. Along the 
"concordance" direction, the flow of the Willick ( 1999) sample is 
820 ± 410 km s~\ 

In contrast, iDale et all J 19991) found a bulk flow of only 
75 ± 92 km s _1 towards I = 289°, b = 25° (with large errors in 
the direction due to the small amplitude of the bulk flow) from TF 
distances in 64 clusters with a similar depth as the SMAC sample 
(but with different spatial sampling). The amplitude is consistent 
with zero. If we fix the direction to be Z = 300°, b = 10°, our fits 
yield an amplitude of 120 ± 120 km s . This is consistent with 
the SMAC flow in this direction. 



6.1.3 EFAR 

IColless et alJ 1200 ll) studied 50 clusters in two distinct regions of 
the sky (Perseus-Pisces and Hercules-Corona-Borealis) via the FP 
method. They found no significant bulk flow for the EFAR sample. 

The geometry of the EFAR survey introduces several sub- 
tle problems. First, there is substantial covariance between the 
monopole (FP zero-point) and dipole terms, with a correlation co- 
efficient of —0.77 between the monopole term and the Galactic 
y-component of the EFAR bulk flow. Since the SMAC bulk flow 
is essentially along the y direction, this is clearly an important is- 
sue. Colless et al. fixed the zero-point of the EFAR FP relation by 
assuming that a subset of their clusters were at rest. However, for 
an anisotropic sample such as EFAR, this choice presupposes that 
there is no bulk flow. They studied the effect of this choice and 
concluded that it was smaller than their random errors. While this 
is correct, it is not negligible and should be included in the total 
errors. Second, irrespective of the monopole covariance, the effects 
of sampling are particularly severe for the EFAR sample, which 
covers two specific regions on the sky. 

We have fit the EFAR peculiar velocities to flow models in- 
cluding both a bulk flow and free zero-point. Note that, because 
we are fitting a flow m odel to distances a posteriori (Method I of 
IStrauss & W illick 1995), we apply homogeneous Malmquist bias 
corrections to the EFAR cluster distances. We also allow for dis- 
tance errors due to errors in the extinctions, in the same way as for 
the SMAC data. Note, however, that the EFAR errors used in this 
analysis do not include the unpublished systematic errors due to 
matching the velocity dispersions systems. We expect these to be 
at a level at least that found for the SMAC sample. 

A fit of the EFAR data yields a bulk flow of 629 ± 38 1 km s " 1 
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towards I = 53° ± 41°, b = 6° ± 25°. If we force EFAR to have 
the same bulk flow as the SMAC sample, we find A^ 2 = 11 for 
an increase of three degrees of freedom. Thus the best fit value 
of the SMAC bulk flow is rejected at the 98.9% CL. This does 
not mean that EFAR results are inconsistent with the SMAC re- 
sults. In the direction I = 300°, b = 10°, the EFAR amplitude 
is 120 ± 310 km s , which is consistent with the SMAC value 
400 ± 120 kms -1 . 



6.1.4 SNla 

We have analyzed 98 SNe within 150 /i _1 Mpc from the compi- 
lation of iTonrv et al] {2003). Details of this analysis will be pre- 
sented elsewhere. In summary, we find that the sample as a whole 
has a bulk flow of 410 ± 75 km s^ 1 towards I = 286°, b = -12°. 
However the effective depth of this sample is only 35/i _1 Mpc. 
To test the coherence of the flow, we have split the sample into 
two subsamples: SNIa-In (0 < r < 60fr _1 Mpc) and SNIa-Out 
(60 < r < 150/jT 1 Mpc). The SNIa-In sample has a bulk flow of 
376 ± 81 km s _1 toward I = 285°, b = -14°, which is not sur- 
prising since this subsample spans a distance range where flows 
are known to be high. The SNIa-Out sample has a bulk flow of 
775±204km s" 1 toward I = 299° ± 17°, b = 3° ± 13°. The bulk 
flow of the outer sample is significantly different from zero at the 
99.8% CL. The direction and amplitude of the SNIA-Out subsam- 
ple are in good agreement with the SMAC bulk flow. 



6.1.5 Comparison of large-scale surveys 

Fig. [5] compares the error-ellipsoids in the V x -V y plane for the 
SMAC, SNIa-Out and Dale et al samples. Note that these errors are 
due to peculiar velocity errors only. They do not include the effects 
of sparse sampling which would tend to further increase the errors. 
These three surveys are the ones which, with the exceptionof LP, 
are most discrepant given their bulk flows and errors. Despite these 
differences there is a "concordance" region in which the flow is 
250-450 km s" 1 . The bulk flows of other surveys mentioned in this 
section cover this concordance region and do not further constrain 
it. 



6.2 Field surveys within approximately dOh Mpc 

The peculiar velocity field has also been mapped at higher density 
within smaller volumes closer to the LG. For these samples, it is 
possible to obtain a bulk flow which is less contaminated by sam- 
pling effects. 

(i) iGiovanelli et"al] (l998) studied a large sample of TF field 
galaxies. They found a bulk flow of 200 ± 65 km s -1 towards 
I = 295°, b = 25°for the volume extending to 6500 km s" 1 . 

(i i) Within 50 ft -1 M pc, the bulk flow of the Mark III sam- 
ple <WillicketalJll997t) is 305 ± 110 kms"' towa rd I = 313°, 
b = 29° when volume- weighted JPekel et af] 1 19991) via the PO- 
TENT method. The corresponding result for the SFI su r vey is 
310 ± 80 kms" 1 toward I = 299°, b = 29°. lDekel etall 119991) 
argue that this flow is generated from beyond 50 ft _1 Mpc. 

(iii) The Shellflow project l Courte au et al J 120001) studied TF 
galaxies in the range 4500 < cz < 7000 km s _1 and found V = 
70±7o° km s" 1 , with a 95% upper limit of V < 300 km s" 1 . Due 
to the small amplitude of the flow in comparison with the errors, the 
direction is not well defined: I = 144° ± 140°, b = 50° ± 79°. 
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Figure 9. Galactic x and y components of the bulk flow and their 68% 
and 95% confidence level error ellipsoids. We have plotted results for the 
large-scale sparse surveys SMAC (solid), Dale et al. (dotted) and SNla with 
60/i Mpc < d < 150/i —1 Mpc (dashed). These surveys all have very 
small components in the z direction. Solutions that lie within the 95% con- 
fidence region of all three surveys are hatched. These three samples which 
differ the most in their flows and hence set the tightest bounds on the concor- 
dance region. The concordance direction I = 300°, b = 10° is indicated 
by the solid diagonal line. The EFAR and Willick results are not plotted as 
their errors are considerably larger. The results of Mould et al. are not plot- 
ted here but are also consistent with this concordance region. In contrast 
to the other surveys, the LP bulk flow vector does have a large z compo- 
nent and so its motion is not well-represented by plotting ellipses in the x-y 
plane. 



(iv) The ENEAR sample of early-type galaxies within cz < 
7000 kms -1 has a bulk flow of V = 220 ± 60 (random) ± 
50 (systematic) km s" 1 toward I = 304° ± 16°, b = 25° ± 11° 
Ida Costa et alfcOOCh. 

tv) ITonrv et al]<2000l) did not measure a bulk flow of their SBF 
sample. Instead, they fitted a combination of internal flows (from 
GA and Virgo) plus external dipole and quadrupole terms. The ex- 
ternal dipole is ~ 150 kms~ toward I — 294°, 6 = 67°. The 
errors on this ampitude are difficult to estimate because of the 
large degeneracies with the GA infall, but are approximately 100 
km s" 1 . 

In summary, with the exception of Shellflow, the field surveys 
find a consistent non-zero flow in a similar direction on the sky. 
Note, however, that the three TF samples (Mark III, SFI, Shellflow) 
are not independent. A conservative estimate of the bulk flow 
within 60 /i _1 Mpc from field TF and FP data is ~ 200±75km s" 1 . 

6.3 Towards a consistent model for large-scale flows 

The bulk flows of all surveys discussed above, projected into the di- 
rection I = 300°, b — 10° are shown in Fie. llQI Surveys are plotted 
as a function of their estimated depth. The lines indicate IRAS PSCz 
predictions for the bulk flow along this direction for different values 
of Pi and V C xt ■ An eyeball fit reveals that V C xt — is excluded by 



Streaming motions of galaxy clusters — V. The peculiar velocity field 15 



XI 500 ; 




5000 



10000 
R (km/s) 



15000 



Figure 10. The bulk flow projected along the direction I = 300°, b = 10° 
as a function of depth, for the LG and va rious surveys: SNIa with r < 
mh- 1 Mpc fSNIa-Irri. lPekel et all tl999l) analysis of the Mark III sam- 
ple (M3-P), SF1, ENEAR, Shellflow, SMA C, distan t (60ft.- 1 Mpc < r < 
150ft- 1 Mpc) SNe (SNIa-out), EFAR and lWillickl <1999h . Error bars are 
lcr and do not include the effects of sparse sampling. Surveys are placed 
approximately at the depth corresponding to the top-hat radius of an ideal- 
ized volume-weighted survey with the same effective depth. This estimate 
is very crude: horizontal error bars on the depth are ~ 10ft- 1 Mpc for the 
nearby samples and ~ 20ft - 1 Mpc for the distant ones. The solid curve 
shows the bulk flow predicted from the PSCz survey if f3\ = 0.5, whereas 
the dotted line shows the predicted bulk flow for /3j = 0.7. The dashed 
lines show the PSCz bulk flow with V cx t of 100, 200 and 300 km s -1 . 



the distant samples, in particular SMAC and the SNIa-Out sample. 
Similarly, a bulk flow with V = 600 km s _1 is also excluded, in 
particular by Dale et ali i 1999h and Shellflow. 

A reasonable fit to all samples is possible if (3i — 0.5, pro- 
ducing a bulk flow component of ~ 100 km s _1 in this direction at 
R = lOO/i" 1 Mpc from sources in the PSCz survey, plus an addi- 
tional Vcxt component from sources not included in PSCz. Most of 
the field peculiar surveys are consistent V cx t of approximately 125 
km s _ , with the exception of Shellflow which prefers V ex t ~ 0. 
The deeper surveys prefer a higher value V cx t ~ 200 kms -1 . 
It is difficult, however, to quantify the effects of sparse sampling, 
particularly on the latter surveys. Inspection of Fig. 1101 suggests 
that V cx t = 125 km s _1 is consistent with almost all of the pecu- 
liar velocity surveys. There appear to be two exceptions: Shellflow, 
which lies too low, and the SNIa-Out sample, which is too high. 
Since both of these samples have unusual geometry — Shellflow 
because it is a shell, and the SNIa-Out sample because it is sparse 
— it is likely that sampling errors will contribute significantly to 
the error bars. 

By adding V cxt = 125 km s" 1 to the 100 km s" 1 
bution from the PSCz density field, one finds that 225 km s 
or a third of the LG's motion arises from sources at large scales 
(7? > 100 /i~ 1 Mpc). Given the systematic uncertainties and effects 
of sampling, we estimate the error on this value is approximately 
100 km s" 1 . 



contn- 

l 



7 CONCLUSIONS 

We have examined in detail the 687 ± 203 km s^ 1 bulk flow of 
the SMAC sample and find that it is robust to systematic errors. We 
have shown that most of the bulk flow is not generated by nearby 
sources such as the Great Attractor but rather arises from structures 
at depths greater than 100 ft _1 Mpc. The Shapley Concentration is 
identified as one likely source of the large scale flow but is unlikely 
to be responsible for all of the SMAC flow. 

When we compare the SMAC motion to the predictions of the 
PSCz survey we find that = 0.39 ± 0.17, consistent with the 
"concordance" value /3 ~ 0.5. However, the IRAS PSCz survey 
can only explain 50 ± 20% of the amplitude of the SMAC flow. 
This suggests that there are gravitational sources not well-mapped 
by the PSCz survey. Evidence from other redshift surveys suggests 
that these sources may be located in the ZOA or in superclusters 
such as Shapley which are undersampled by IRAS. 

Finally, we have shown that the SMAC survey is not inconsis- 
tent with oth er sparsely-sampled survey s of the large scale velocity 
field, such aslDale et all h999h : lwiiiickl<l 9991 EFAR) and the SNe 
sample of iTonrv et al. 1 2003). Taken together, all surveys suggest 
a large scale flow approximately 225 kms -1 towards / = 300°, 
b — 10°. Further analyses of the effects of sparse sampling, and 
detailed comparisons with the predictions of redshift surveys such 
as PSCz are needed to more accuractely quantify this result. 

The NO AO Fundamental Plane Survey Ismith et alj |2004 
will measure the peculiar velocities of 100 X-ray selected clus- 
ters. The total number of FP distances will be ~ 4000, about six 
times the number in SMAC. When results from the NFPS are com- 
pared to, for exam ple, predictions of the 2MASS redshift survey 
iHuchraetail E003), we expect to identify the sources the LG's 
motion. 
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